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15 Abstract 

16 We examine the ozone production from boreal forest fires based on a case study of 

17 wildfires in Alaska and Canada in summer 2004. The model simulations were performed 

18 with the chemistry transport model, MOZART-4, and were evaluated by comparison with 

19 a comprehensive set of aircraft measurements. In the analysis we use measurements and 

20 model simulations of carbon monoxide (CO) and ozone (O3) at the PICO-NARE station 

21 located in the Azores within the pathway of North American outflow. The modeled 

22 mixing ratios were used to test the robustness of the enhancement ratio AO 3 /ACO 

23 (defined as the excess O3 mixing ratio normalized by the increase in CO) and the 

24 feasibility for using this ratio in estimating the O3 production from the wildfires. Modeled 
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1 and observed enhancement ratios are about 0.25 ppbv/ppbv which is in the range of 

2 values found in the literature, and results in a global net O 3 production of 12.9±2 Tg O 3 

3 during summer 2004. This matches the net O 3 production calculated in the model for a 

4 region extending from Alaska to the East Atlantic (9-11 Tg O3) indicating that 

5 observations at PICO-NARE representing photochemically well-aged plumes provide a 

6 good measure of the O3 production of North American boreal fires. However, net 

7 chemical loss of fire related O 3 dominates in regions far downwind from the fires (e.g. 

8 Europe and Asia) resulting in a global net O3 production of 6 Tg O3 during the same time 

9 period. On average, the fires increased the O3 burden (surface-300 mbar) over Alaska 
10 and Canada during summer 2004 by about 7-9%, and over Europe by about 2-3%. 
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1 


1 Introduction 


2 Ozone (O3) plays a central role in tropospheric chemistry as a primary source of 

3 hydroxyl radicals and, by being toxic in nature, has negative impacts on human and plant 

4 health. It is also estimated to be the third most important anthropogenic greenhouse gas 

5 [Ramaswamy et al., 2001], Anthropogenic sources and biomass burning release O 3 

6 precursors including carbon monoxide (CO), nitrogen oxides (NO x ) and volatile organic 

7 compounds (VOCs) into the atmosphere. Photochemical reaction of CO and VOCs with 

8 the hydroxyl radical in the presence of NO x and sunlight results in the production of O3. 

9 The production of tropospheric ozone in the Northern mid-latitudes is largely 

10 impacted by anthropogenic sources [ Chameides and Tan, 1981; Levy et al., 1985]. 

1 1 Significant ozone enhancements have been observed in individual plumes of boreal forest 

12 fires [Wofsy et al., 1992; Goode et al., 2000; Forster et al., 2001; McKeen et al., 2002| 

13 Jaffe et al., 2004; Honrath et al., 2004; Lapina et al., 2006], and measurements in 

14 combination with chemical transport simulations have been used in various studies to 

15 estimate the amount of ozone produced from boreal fires [Mauzerall et al., 1996; 

16 McKeen et al., 2002], However, the large-scale impacts of high latitude biomass burning 

17 on the hemispheric tropospheric ozone budget are poorly quantified. 

18 In here we apply a chemical transport model to evaluate various techniques used for 

19 estimating the ozone production from a specific source, and include various model tracers 

20 to gain a detailed insight into the limitations of these methods. We combine model 

21 analysis with observations of CO and O3 to quantify contributions of boreal fires to 

22 Northern Hemispheric CO and O3 burdens, a topic not very well explored so far. Our 
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1 analysis focuses on fires in Alaska and Canada in summer 2004. These fires were the 

2 largest on record for Alaska, and the CO emissions for the North American boreal region 

3 has been estimated as 30±5 Tg for June through August [Pfister et al., 2005]. A total of 

4 about 1 1 million acres were burned in Alaska and Canada during that time period. The 

5 study is supported by a comprehensive set of observations collected during the ICARTT 

6 (International Consortium for Atmospheric Research on Transport and Transformation) 

7 campaign. 

8 CO is a long-lived tracer, and the relationship between mixing ratios of O 3 and CO in 

9 transported regional plumes can be used as an indicator for the magnitude of net O 3 

10 production from selected sources [Parrish et al., 1993], It has been found that the 

1 1 enhancement ratio (AO 3 /ACO), given as the excess O 3 mixing ratio normalized by the 

12 increase in CO concentrations, is typically smaller for boreal forest fires than for tropical 

13 biomass and savannah burning or urban and industrial plumes, due to a lower NO x :CO 

14 emission ratio in boreal forest fires compared to the other sources [e.g. Andreae et al., 

15 1994; Wofsy et al., 1992], AO 3 /ACO of fire plumes is also expected to change with plume 

16 age. For example, Yokelson et al. [2003] found an increase from 0.09 ppbv/ppbv in fresh 

17 tropical biomass burning plumes to 0.22 ppbv/ppbv for plumes 2-4 days old. Thus, O 3 

1 8 production downwind from the source region must be accounted for. 

19 The structure of this paper is the following. After the Introduction we discuss the 

20 model simulations and model evaluation in Sections 2 and 3, respectively. In Section 4 

21 we describe CO and O3 in-situ measurements taken at the PICO-NARE station located in 

22 the Azores. These observations were used in combination with model simulations to 

23 investigate the O 3 production from the fires in Alaska and Canada in summer 2004. 
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1 Section 5 discusses and evaluates different techniques for calculating the enhancement 

2 ratio and analyzes the O 3 production due to emissions from the fires. The analysis is 

3 supported by incorporating fire tracers for CO and O3 into the model and by performing 

4 model simulations with and without fire emissions. Finally, we investigate the 

5 contributions these fires had on the Northern Hemispheric and regional budgets of CO 

6 and O 3 . Section 6 summarizes our findings. 
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2 Model Simulation 


The Model for OZone And Related chemical Tracers (MOZART) chemistry transport 
model has been developed at the National Center for Atmospheric Research, the 
Geophysical Fluid Dynamics Laboratory and the Max-Planck Institute for Meteorology. 
In this study we are using Version 4 [any other reference for Mozart with aerosols?, 
Emmons et al., Sensitivity of chemical budgets to meteorology in MOZART-4, in 
preparation]. Modifications from Version 2 published in Horowitz et al. [2003] include, 
amongst others, a more complete description of anthropogenic hydrocarbon chemistry, 
the inclusion of tropospheric aerosols, and on-line calculations of dry deposition, FLO, 
and biogenic emissions. 

We run the model at a horizontal resolution of ~2.8 degrees by 2.8 degrees. The 
meteorological fields for 2004 for driving MOZART were taken from the National 
Centers for Environmental Prediction (NCEP) National Center for Atmospheric Research 
Re- Analysis [ Kistler et al., 2001] and were interpolated from a 6-hour time resolution to 
the 20-minute time steps of the simulations. The vertical resolution of the meteorological 
fields and hence the model consists of 28 hybrid levels ranging from the surface up to 2 
hPa. 

Biofuel and fossil fuel emissions used in this study were taken from the European 
Union project POET (Precursors of Ozone and their Effects in the Troposphere) [Granier 
et al., 2004], Over the continental US, the anthropogenic emissions are based on the U.S. 
EPA NEI-99 inventory (National Emissions Inventory, base year 1999, version 3) [EPA, 
2004], For the Alaska and Canada region, the biomass burning emissions for CO for 2004 
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1 


were taken from an inverse modeling study [ Pfister et al., 2005], and emissions for NO : 


2 and VOCs were deduced from this inventory by applying emission factors based on 

3 Andreae and Merlet [2001], At the time these simulations were run, an emissions 

4 inventory for the year 2004 for biomass burning sources outside North America was not 

5 available. A comparison of CO data from the Measurements Of Pollution In The 

6 Troposphere (MOPITT) remote sensing instrument for 2000 - 2004 showed that the 

7 global biomass burning activity in summer 2004 was similar to 2002 and for this reason 

8 we used a 2002 biomass burning inventory based on ATSR fire counts [Granier et al., 

9 2004], 

10 Our model simulations cover the months from June through August 2004 with a spin- 

11 up phase beginning in August 2003. We performed three different simulations. Two of 

12 these include emissions from the Alaskan and Canadian wildfires, and are abbreviated as 

13 BB in the following. In one of these simulations ( BBsrJ ), the wildfire emissions were 

14 released at the lowest model layer and distributed in the boundary layer by the model 

15 boundary layer scheme. Studies [Fromm et al., 2005; Damoah et al., 2006] have shown 

16 that that fire induced convective clouds might transport fire emissions rapidly to higher 

17 altitude. To test the sensitivity of our model to the emissions injection height, we 

18 performed another simulation (BBvert), where the emissions were distributed evenly with 

19 regard to number density between the surface and 9 km altitude. The 9 km altitude 

20 represents an upper limit for the injection height based on estimates derived from the 

2 1 Multi- Angle Imaging Spectro-Radiometer (MISR) and the Moderate Resolution Imaging 

22 Spectro-Radiometer (MODIS) [Averill et al., 2005]. In a third simulation (noBB) used as 

23 a reference, the emissions of the wildfires in Alaska and Canada were set to zero. 
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1 We included two fire tracers in the model. Tracers are emitted or produced from a 

2 specific source, but undergo the same transport, chemistry, and physical processes as the 

3 standard species. The first tracer incorporated into the simulations is a CO fire tracer 

4 (COJ), i.e. CO released from the Alaska/Canada wildfires. For the second tracer we 

5 tagged the O 3 production resulting from hydrocarbon or CO oxidation in association with 

6 the emissions of NO x from the fires. We refer to this tracer as Oj NOx in the following. The 

7 tagging technique for Os x takes into account the re-cycling of NO x from reservoirs such 

8 as PAN by applying tags to all nitrogen-containing species. Although there are some 

9 minor pathways to create O3 without the presence of NO x , the accuracy of the tagging 

10 technique has been estimated as better than 95% on a monthly basis [ Lamarque et al., 

1 1 2005], The statistical analysis we perform using Os NOx in this study is expected to give a 

12 comparable accuracy. 
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1 


3 Model Evaluation 


2 MOZART-4 simulations ( BBsrf and BBvert ) have been evaluated by comparison with 

3 aircraft measurements taken in the framework of the ICARTT campaign during summer 

4 2004. Table 1 includes a list of the platforms, instrumentation, and corresponding 

5 references. Measurements with the NASA DC-8 covered large parts of North America 

6 and the Northern Atlantic, the NOAA-P3 flights were focused on the Eastern United 

7 States, the British BAE146 performed most flights over the Atlantic, and the German 

8 Falcon covered mostly Europe. The Measurements of OZone aboard Airbus In-service 

9 airCraft (MOZAIC) data set has global coverage, but we are including only 

10 measurements over North America, the Northern Atlantic and Europe in accordance with 

1 1 the regions covered by the other aircraft. Most low altitude measurements for this data set 

12 are from take-offs and landings over airports in Europe and the US. For a more detailed 

13 description of all flight patterns we refer to Fehsenfeld et al. [International Consortium 

14 for Atmospheric Research on Transport and Transformation (ICARTT): North America 

15 to Europe: Overview of the 2004 summer field study, submitted to J. Geophys. Res., 

16 2006], 

17 For the comparison of the model with the observations, 2-hour average model data 

18 have been linearly interpolated to the time and location of the aircraft data. The time 

19 resolution of the observations is 1 minute. Statistics have been calculated for the 

20 individual aircraft data sets and binned onto a 2-km wide altitude grid. The results for 

21 BBvert, BBsrf, and for comparison, noBB , are shown in Figure 1 for CO and in Figure 2 

22 for O3, respectively. The agreement with observations for CO as well as O3 is generally 
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1 better in the simulations with fires than in the simulation without fires. CO and O 3 mixing 

2 ratios are clearly different between the BB and noBB simulations for all platforms, also 

3 for the UK BAE146 and DLR Falcon. This indicates that plumes from the Alaska and 

4 Canada fires reached all the way to Europe. We list statistics for modeled CO and O3 

5 concentrations for the different platforms in Table 2. The data included in the statistics 

6 corresponds to those shown in Figures 1 and 2, but in here we consider the entire altitude 

7 range from the surface to 8 km. The T-test significance levels for simulations with and 

8 without fire emissions are above 99% for all platforms indicating the samples have 

9 significantly different means. It is interesting to note that the T-test statisctics comparing 

10 BBvert and BBsrf show a somewhat higher significance level for O3 compared to CO. 

1 1 This suggests the modeled O 3 production from the fires is slightly more sensitive to the 

12 injection height than the concentration fields of CO. 

13 For most altitudes and platforms, the CO bias between model and observations is less 

14 than 10% for both BBsrf and BBvert (see Figure 1). The mean bias as well as the 

15 correlation improved upon adding fire emissions into the model, with the only exception 

16 being the highest altitude bin for the NOAA P3 and DLR Falcon data set. We believe this 

17 can partly be explained by the small data sample in these bins, and, associated with that, 

18 the comparison is more strongly impacted by single events. This might also contribute to 

19 the large bias in the 4-6 km bin for the DLR Falcon. The comparison for O3 (Figure 2) 

20 shows an agreement of better than 10 % for all platforms and altitude bins when fire 

21 emissions are included in the simulations. No clear conclusion can be drawn from the 

22 evaluation if either BBsrf or BBvert lead to better general agreement. The comparison for 

23 individual fire plumes give a better agreement for either the one or the other simulation. 
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1 This reflects the combination of crown, smoldering, and peat burning of the Alaska fires. 

2 The two cases we ran do not represent the full complexity of fire behavior, but are 

3 probably better regarded as sensitivity tests to the vertical distribution of emissions. 
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1 


4 Impact of Biomass Burning on CO and 0 3 at PICO-NARE 


2 For the analysis of the O3 production from the wildfires in Alaska and Canada we 

3 made use of in-situ measurements at the PICO-NARE station. The station is located on 

4 the summit caldera of Pico mountain on Pico Island in the Azores, Portugal (2225 m 

5 above sea level, 38.47N, 28.40W) and is well suited for studying North American 

6 pollution outflow. Air masses at this location typically arrive from North America, but 

7 frequently originate from high latitude regions such as Alaska and Siberia [Honrath et 

8 al, 2004], often with enhancements in CO and O3 that have been attributed to boreal fire 

9 impacts [Lapina et al, 2006], The advantage of measurements at PICO-NARE compared 

10 with other locations on the continent is its remote location allowing the sampling of 

1 1 chemically well processed air masses. 

12 CO at PICO-NARE was measured using a non-dispersive infrared absorption 

13 instrument (Thermo Environmental, Inc., Model 48C-TL) modified as described by 

14 Parrish et al. [1994]. O3 was measured with a commercial ultraviolet absorption 

15 instrument (Thermo Environmental Instruments Inc., Franklin, Massachusetts, Model 

16 49C). Data are available as 1 -minute averages. For the 2-hour averages used in this study 

17 the precision for CO is estimated to be better than 9 ppbv and for O3 better than 1 ppbv. 

18 For a description of the station and the measurement techniques we refer to Honrath et al. 

19 [2004] and Owen et al. [2006]. 

20 For comparing the model simulations with the observations at PICO-NARE, the 

2 1 observations have been averaged in time to match the 2-hour window of the simulations, 

22 and the model data have then been linearly interpolated to the location and pressure level 
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1 


of the observations. Due to the coarse model resolution, the simulations do not resolve 


2 upslope events occurring at the mountain site. Observations potentially affected by 

3 upslope flow were small during summer 2004. We identified these periods as described 

4 by Kleissl et al. [2006] [The occurrence of upslope flows at the Pico mountain-top 

5 observatory: a case study of orographic flows on a small, volcanic island, submitted to J. 

6 Geophys. Res.] and omitted them from the analysis. 

7 Figure 3 shows the time series for modeled and measured CO and O3 at P1CO-NARE. 

8 In addition to results from the BB model simulations, we also include results from the 

9 noBB simulation to emphasize the impact of the Alaskan/Canadian wildfires. For clarity, 

10 we reduced the temporal resolution in the graphs to daily average values, however, our 

1 1 analysis refers to the 2-hour average values. A detailed analysis of P1CO-NARE 30- 

12 minute observations in the 2004 fire season is provided by Val Martin et al. [Significant 

13 Enhancement of Nitrogen Oxides, Black Carbon, and Ozone in the North Atlantic Lower 

14 Free Troposphere Resulting from North American Boreal Wildfires, submitted to J. 

15 Geophys. Res. (hereinafter Val Martin et al., submitted)]. As can be seen in Figure 3, the 

16 model closely matches the observed temporal variability, and also captures the 

17 magnitudes fairly well. The mean bias between modeled ( BBsrf) and observed mixing 

18 ratios is -3±16 ppbv for CO and 6± 12 ppbv for O3. The corresponding biases for the 

19 simulation without fire emissions are — 12± 1 9 ppbv for CO and 4± 12 ppbv for O3, 

20 respectively. Daily CO fire tracers in the model estimate a typical transport time on the 

21 order of 1-2 weeks for biomass burning plumes reaching P1CO-NARE. 

22 The time series for CO from the simulations BB and noBB indicate a clear impact from 

23 the fires at P1CO-NARE throughout the months of July and August with the most 
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1 extreme period from the end of July to mid-August. This is a combination of the most 

2 intense fire activity occurring in mid to end of July [ Pfister et al., 2005] and the transport 

3 time of 1-2 weeks. The period from July 22-24 had the highest half-hour average CO 

4 levels yet recorded at the PICO-NARE station [ Val Martin et al, submitted]. The impact 

5 of the fires on the O3 concentrations is less pronounced, but differences between BB and 

6 noBB of up to 10 ppbv are evident during some of the intense episodes. The difference in 

7 CO and O3 concentrations by subtracting noBB from BBsrf gives an average 

8 enhancement due to the fires of 8 ppbv (8%) for CO and 2 ppbv (4%) for O3. 

9 The correlation between the CO mixing ratios from the simulation BBsrf and the 

10 observations is r=0.64 compared to a correlation of r=0.48 between noBB and the 

1 1 observations. During times of intense biomass burning impact, the noBB run actually 

12 shows slight enhancement in the CO concentrations as well, indicating that these outflow 

13 events transported pollution from the fires together with elevated pollution from likely 

14 North American anthropogenic sources. 

15 The correlation between the measured and modeled O3 is r=0.51 for the simulation 

16 noBB and increases only slightly for the BB runs (r=0.54), an indication of the less 

17 pronounced or more complex [ Val Martin et al., submitted] effect of the fires on the O3 

1 8 burden compared to the CO burden. While the model is picking up the higher O3 values 

19 moderately well, neither the BB nor the noBB simulations capture the low end of the 

20 observed O3 concentrations. This is likely due to the positive O3 bias in the model over 

21 the US, leading to an overestimate of North American outflow of O3 and to an 

22 overestimate in the modeled O3 mixing ratio of maritime background air. In agreement 

23 with the conclusions drawn from the model evaluation in Section 3 it cannot be stated if 
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1 either BBsrf or BBvert results in better agreement with the observations. Unless otherwise 

2 mentioned, we will focus in the following analysis on results for BBsrf. 
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5 Ozone Production from Boreal Fires 


2 Assuming a linear relationship between CO emissions and net O3 production, the 

3 relationship between tropospheric CO and O3 concentrations might be used as an estimate 

4 for the net O3 production in regional plumes [Parrish et al., 1993; Forster et al., 2001], 

5 The enhancement ratio AO3/ACO is defined as the difference between the O3 

6 concentrations in a polluted air mass from that of background air, normalized by the 

7 excess mixing ratio of CO. In the case of biomass burning plumes, the background 

8 defines concentrations of CO and O3 not linked to the fire emissions. 

9 The change in the Northern Hemispheric net O3 production rate due to the fires, 

10 calculated by differencing net O3 production rates in the BB and the noBB simulation, is 

1 1 estimated as 6 Tg O3 for June through August. This is contribution of 3% to the Northern 

12 Hemispheric net ozone production. Normalizing by the total CO emissions for this time 

13 period (30±5 Tg CO) yields a global average enhancement ratio for the fires in Alaska 

14 and Canada of 0.12 ppbv/ppbv. In the following Sections we examine the feasibility of 

15 estimating the O3 production of North American boreal fires by using enhancement ratios 

16 based on observed and modeled mixing ratios of CO and O3 at PICO-NARE. 

17 5.1 Determining the Enhancement Ratio from CO and O3 Observations 

18 There are two common ways to calculate AO3/ACO. The first, in the following 

19 referred to as the “scatter technique”, determines the enhancement ratio from the slope of 

20 the linear fit of O3 versus CO mixing ratios [Parrish et al., 1993]. The second, termed as 

21 “enhancement technique”, infers “background” concentrations of CO and O3 from air 
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masses not affected by the considered pollution source (in our case the wildfires in 


2 Alaska and Canada) and calculates the corresponding excess mixing ratios by subtracting 

3 background concentrations from total concentrations [Mauzerall et al., 1998], Both 

4 methods are discussed in the following. 

5 5.1.1 Scatter Technique 

6 This technique has been applied in Figure 4 showing CO-O 3 scatter plots for observed 

7 and modeled concentrations at PICO-NARE. The data were grouped into air masses with 

8 varying biomass burning impact by using the magnitude of the observed CO mixing ratio 

9 as the threshold. Studies by Honrath et al. [2004] and Lapina et al. [2006] show that 

10 periods of extreme summertime CO concentrations frequently coincide with airflow from 

1 1 Northern latitudes transporting pollution from wildfires in Siberia, Alaska and Canada to 

12 PICO-NARE. To allow comparison with the observations, we applied two methods to the 

13 model data: 1) we used a threshold derived from the simulated CO, and 2) we used the 

14 relative contribution of the fire tracer COf. The second method gives in some sense the 

15 true solution, as it is not impacted by sources of high CO other than the wildfires. The 

16 thresholds, specified in the graphs, were chosen in a way that the number of data points 

17 was roughly equal in the different subsets. The coarse spatial resolution in the model 

18 results in a more pronounced dilution of biomass burning plumes and this explains the 

19 smaller threshold in total CO applied to the subset of high intense plumes in the model 

20 compared to the observations. 

21 The fitting technique applied is a reduced major axis reduction (RMA). It uses the 

22 geometric mean of the slopes of the standard linear regression of y versus x and of x 

23 versus y [Draper and Smith, 1998] thus taking into account the variability in both 
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1 abscissa and ordinate. The mean slope and the corresponding standard deviation are 

2 specified in Figure 4. 

3 The scatter in the data and the uncertainties in the slopes are explained by the mixture 

4 of different air masses, variability in the background CO and O 3 levels, and differing 

5 pathways and photochemical ages in the sampled plumes. The slopes are also somewhat 

6 dependent on the fitting technique applied. The calculated AO3/ACO is independent of air 

7 mass age if the tracers used have similar lifetimes or have lifetimes much longer than the 

8 transport time, but this assumption is not entirely fulfilled in the case of CO and O3. 

9 However, in a statistical sense, a clear distinction between the different types of air 

10 masses is evident with the smallest O3 enhancements per unit CO for the most intense 

1 1 plumes in both model and data. The decrease in the slopes with increasing biomass 

12 burning impact is seen in both the model data (Figure 4b) and the observations (Figure 

13 4a). However, the model does not capture the measured AO 3 /ACO due to the impact of 

14 mixing in the model. As a result, the calculated slope has a rather large value with, at the 

15 same time, a high uncertainty. Applying the ratio COf/CO as a selection criteria instead 

16 (Figure 4c) we achieved a clearly stricter separation of the fire impact for the modeled 

17 plumes. The slope calculated for the most biomass burning impacted subset is then 

18 0.44±0.08 with a correlation of r=0.54 instead of a slope of 0.77±0.57 and a correlation of 

19 r=0.28 when applying total CO mixing ratios as threshold. This slope is still higher than 

20 the observed one and this is likely due to more strongly diluted plumes in the model 

21 resulting from the coarse spatial resolution. As will be shown in Section 5.2, the selection 

22 of more intense model fire plumes reduces the estimated slope. 


Pfister_INTEX_revision.doc Version from 9/9/06 


Page 18/45 



1 The generally higher slopes for air masses least impacted by biomass burning sources 

2 (defined as “non-fire plumes” here and characterized strongly by anthropogenic pollution 

3 sources) compared to those more strongly impacted by the wildfires is consistent with 

4 earlier studies [Wofsy et al., 1992; Andreae et al., 1994, McKeen et al., 2002] and have 

5 been explained by a lower NO x :CO emission ratio in boreal fires relative to urban and 

6 industrial sources In Figure 4c we include results from a simulation where the NO x :CO 

7 emission ratio for the wildfire emissions in the model was increased by a factor of 10 to 

8 match the emission ratio of anthropogenic sources. The pronounced difference between 

9 air masses with weak and with strong biomass burning impact is diminished in this case 

10 supporting the hypothesis that the difference between O 3 production of anthropogenic 

1 1 and biomass burning plumes can largely be explained by a difference in NO x :CO 

12 emission ratios. 

13 It is evident from Figure 4 that the limitations used in the selection of biomass burning 

14 impacted air masses and the degree of mixing in the considered air masses have an effect 

15 on the calculation of AO 3 /ACO. The observed and modeled correlations also indicate that 

16 the enhancement ratio decreases the more strict the selection criteria applied are, and this 

17 dependence will be looked into more closely in Section 5.2. 

18 5.1.2 Enhancement Technique 

19 The enhancement technique requires knowledge of the background concentrations, 

20 and large uncertainties might be introduced if the background and its variability are not 

21 well known. For modeling studies, however, the variability in background concentrations 

22 can be determined accurately from a simulation where the considered emission source is 

23 omitted. 
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1 Figure 5 shows observed and modeled excess mixing ratios of O 3 versus CO for the air 

2 mass separation applied earlier (Section 5.1.1). For Figures 5a and 5b we estimated the 

3 background by averaging CO and O3 mixing ratios over the subset of non-fire impacted 

4 air flow, which explains the existence of negative excess mixing ratios in this subset of 

5 data. For the observations the derived background values are 82 ppbv for CO and 37 

6 ppbv for O 3 ; for the model data we obtain 81 ppbv for CO and 45 ppbv for O 3 . The 

7 derived slopes for the air masses impacted by biomass burning are nearly identical to the 

8 values derived from the scatter technique showing that the background value we assumed 

9 in the enhancement technique is similar to the one implied in the scatter technique. The 

10 difference in the slopes for observed and modeled CO-O 3 relationships (Figures 5a and 

1 1 5b) suggests this method is also sensitive to the amount of mixing. 

12 Figure 5c shows modeled slopes when the background is estimated from the 

13 simulation without fire emissions. In some sense, this is the result one would obtain 

14 assuming perfect data, i.e. if the excess mixing ratios of CO and O3 are precisely known. 

15 For each data point in the BBsrf simulation a corresponding background value is derived 

16 from the noBB run, thus the temporal variability in the background is accounted for. The 

17 standard deviation of the background mixing ratios as determined from the noBB 

18 simulation is on the order of 5 ppbv for the biomass burning plumes. Compared to the 

19 scatter technique and the enhancement technique with a constant background, the excess 

20 mixing ratios for CO and O3 derived with this technique are more strongly correlated and 

21 the calculated slopes are much less sensitive to the selection of the air mass. AO3/ACO for 

22 the biomass burning plumes at PICO-NARE derived from this technique is 0.28 

23 ppbv/ppbv when the air masses are filtered for the points most impacted by biomass 
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1 burning, and 0.23 ppbv/ppbv for all data points. These values are lower than those 

2 calculated with the scattering technique for the model data and match the observed 

3 AO3/ACO. This concludes that the modeled fire plumes are more diluted compared to 

4 observed plumes because of the coarse model resolution, and still carry characteristics of 

5 non- fire related pollution. This impact only cancels out by considering the contribution of 

6 time-varying non-fire related background air. When the sensitivity to model mixing is 

7 removed, the model captures the observed enhancement ratio. 

8 Thus, provided accurate information about the variability in the background levels is 

9 available, the enhancement technique allows a more accurate determination of AO3/ACO 

10 than the scatter technique. The same technique applied to the model simulation in which 

1 1 the NO x fire emissions were increased by a factor of 10 yields slopes on the order of 0.8, 

12 i.e. close to the enhancement ratio estimated for anthropogenic sources. 

13 5.2 Enhancement Ratios and O 3 Production 

14 We used the model data to test the sensitivity of AO3/ACO derived with both the 

15 scatter and the enhancement techniques to the degree to which air masses are impacted by 

16 biomass burning. For the enhancement technique the background was derived from the 

17 noBB simulation. Figure 6 shows AO3/ACO as a function of the magnitude of biomass 

18 burning impact indicated by selecting air masses based on a lower limit of the fraction of 

19 COf to total CO. The enhancement technique shows a weak dependence on the selected 

20 air masses and the variations seen reflect the variability in O3 production in biomass 

21 burning plumes due to changes in O3 chemistry, transport pathways and times, and a 

22 combination of these or more processes. The average AO3/ACO is calculated as 0.25 
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1 


ppbv/ppbv for BBsrf and a slightly higher value, 0.27 ppbv/ppbv, is calculated for 


2 BBvert. This value is close to the enhancement ratio of 0.26 derived from the 

3 observations when intense plumes are selected only (Figure 4). 

4 Our values for AO3/ACO are in the range of values found in the literature. Mauzerall 

5 et al. [1996] calculated enhancement ratios for aged boreal fire plumes on the order of 

6 0.1 ±0.2 ppbv/ppbv. A value of 0.1 ppbv/ppbv was encountered during SOS-95 by 

7 Wotawa and Trainer [2000] and similar values during ABLE-3 for Alaska fires [Jacob et 

8 al., 1992], McKeen et al. [2002] report enhancement ratios of 0.17 ppbv/ppbv. Higher 

9 averaged enhancement ratios are estimated by Bertschi and Jaffe [2005] and Honrath et 

10 al. [2004] for highly aged boreal fire plumes: 0.4 ppbv/ppbv and 0.7 ppbv/ppbv, 

1 1 respectively. 

12 In the case of the scatter technique it is evident that AO3/ACO is higher for the weakly 

13 impacted plumes due to the mixing of the biomass burning impact with the impact of 

14 other pollution sources. As mentioned earlier, during times of intense fire plumes, 

15 increased pollution was also transported to PICO-NARE (Figure 3). The average O3 

16 concentration for the noBB simulation is 44±9 ppbv for COf/CO< 0.01, but 47±5 ppbv 

17 and 47±4 ppbv for C()f/CO>0. 1 and COf/CO> 0.2, respectively. Corresponding CO 

18 concentrations are 77±8 ppbv, 84±6 ppbv, and 82±6 ppbv. When air masses with at least 

19 20% biomass burning impact are selected, the slopes calculated with the two different 

20 techniques approaches a similar range. However, towards stricter limitations the number 

21 of data points is small and the fitting technique is less reliable. 

22 By using AO3/ACO = 0.25 ppbv/ppbv as derived from the enhancement technique, an 

23 approximation for the total O3 produced from the fires can be made [. Parrish et al., 1993; 
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1 Mauzerall et al., 1996], With 30±5 Tg CO emitted by the fires from June through August 

2 as stated by Pfister et al. [2005], an O 3 production of 10.7-15 Tg O 3 is estimated for the 

3 same time period (30±5 Tg CO multiplied by 0.25 ppbv ppbv ' 1 and corrected by the ratio 

4 of O3 (48 g mol' 1 ) to CO (28 g mol' 1 ) molecular weights). 

5 This value is larger than the Northern Hemispheric net chemical production rate of 

6 6 Tg O3 in the model as mentioned earlier. The discrepancy can be explained in that 

7 AO3/ACO at the location of PICO-NARE is not representative for the total net change in 

8 O3, but rather for the net O3 production rate covering the region from the source location 

9 to the Azores. Close to the source region net chemical production of fire-related O 3 

10 dominates, while further downwind from the source (e.g. Europe and Asia) net chemical 

1 1 loss dominates. The modeled net O 3 production calculated over a region representative 

12 for air masses reaching PICO-NARE (stretching from 180W to 20W and from 40N to 

13 70N) is 9 Tg O3 for BBsrf and 1 1 Tg O3 for BBvert, i.e. in the range of the estimate based 

14 on AO 3 /ACO. These results demonstrate that the measurements and model simulations at 

15 the location of PICO-NARE being representative of aged biomass burning plumes indeed 

16 give a good measure of the O 3 production of fires in North America. Differences between 

17 the AO 3 /ACO based estimate and the model calculated production are explained by 

18 uncertainties in the calculation of the slopes and differing pathways and chemical ages 

19 for plumes reaching PICO-NARE. 

20 As mentioned earlier, our model simulations not only include a CO fire tracer, but also 

21 an O3 fire tracer 0/ :(h that tracks the amount of O3 produced from the NO x fire 

22 emissions. One might assume that the global net O3 production for this tracer equals the 

23 amount derived when subtracting results for simulations with and without fire emissions 
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1 as has been done above (6 Tg O3). However, the net O3 production rate for the Northern 

2 Hemisphere calculated from the O3 fire tracer is higher, close to 9 Tg O3. The reasons for 

3 this are due to the non-linearity in O3 chemistry and are explored in the following 

4 Section. 

5 5.3 Changes in the O 3 Chemistry due to Fire Emissions 

6 In Figure 7 we show correlations between the model CO and O3 fire tracers and the 

7 difference in O3 and CO mixing ratios from BBsrf and noBB simulations defined as <70; 

8 and dCO, respectively. As can be seen, the correlation for CO is close to the 1:1 line, but 

9 for O3 the mixing ratios of Of'° x are clearly larger compared to d0 3 indicating that O3 

10 production related to the fires must have been offset by an increased loss of O3. 

11 To explore the mechanisms behind the O3 production from the fire emissions, we 

12 compared O3 concentrations, production and loss terms in the model for the simulations 

13 BBsrf and noBB. The data set has been split into three groups of varying fire impact 

14 determined by the ratio COf/CO. Statistics for the individual subsets are plotted in Figure 

15 8. The maps (Figure 8a) denote the geographical coverage of the selected data with red 

16 indicating a high, and blue a low, concentration of data points. The most intense plumes 

17 are concentrated near the source location, but plumes of high fire impact can also be seen 

18 all the way to Europe. The high intensity plumes are mostly located near the surface 

19 (Figure 8b) as expected when the emissions are released at the lowest model level. With 

20 time, atmospheric transport and convection spread plumes over a larger altitude range. 

21 Figure 8c shows that O3 levels without fire emissions (OpioBB) are mostly below 40 

22 ppbv over the source regions, i.e. the fires occurred in an area of low O3 concentrations. 
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1 The additional NO x from the fires causes a shift in the distribution of O 3 concentrations 

2 towards higher values. The effect is most pronounced over the source region, but a slight 

3 positive shift is also evident for the subset of least impacted plumes. 

4 The increase in O 3 is caused by a strong net production of the fire tracer 0/ vr; ' in the 

5 most impacted plumes with the magnitude decreasing with decreasing plume intensity. 

6 This is shown in Figure 8 d where we illustrate histograms for the C)/' (h net chemical 

7 production. Even though the strongest production takes place close to the source region, 

8 continuing production is also evident in regions further downwind from the source. For 

9 the subset of least impacted plumes (and most aged plumes), there are a significant 

10 number of data points with net chemical loss of C> 3 NOx . 

1 1 Figures 8 e to 8 g show the changes to background O3 levels when fire emissions are 

12 injected into the system. Figure 8 e denotes the distribution for production of O 3110 BB and 

13 of background O3 with fire emissions; the latter defined as O3B. O3B is calculated by 

14 subtracting the O 3 fire tracer from the total O 3 . The production of O3B is less than the 

15 production of 0 3 noBB. Thus, by adding fire emissions to the system, the production of 

16 background O 3 is reduced in the simulations. This is explained by reduced background 

17 levels of peroxy radicals (not shown here). Peroxy radicals play an important role in 

18 ozone production by reacting with NO to form NO 2 which then is photolyzed to give 

19 atomic oxygen necessary in the O3 formation. Changes are most pronounced over the 

20 source region, but differences are also evident in less impacted plumes. 

21 In addition to a reduced production of O3B, the loss of O3B is increased over that of 

22 O 3110 BB, reflected in a reduction in the chemical lifetime (Figure 8 f). The reduced 

23 production and increased loss of background O 3 result in lower concentrations of O3B 
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1 compared to OsnoBB. In Figure 8g we show the corresponding frequency distributions 

2 for Os, OsB and dOs . When compared to Figure 8c we see that over the source region 

3 OsB has on average -80% smaller values compared to OsnoBB and most of the O3 

4 present is in the form of Os NOx , i.e. O3 due to the NO x fire emissions. 

5 The strong reduction of background O3 levels when fire emissions are included 

6 explains why concentrations of Os NOx are larger than the difference in O3 concentrations 

7 between runs with and without fire emissions. In contrast, the CO chemistry has a first- 

8 order linearity, thus the concentrations of the fire tracer COf are close to the difference of 

9 CO concentrations simulated with and without fire emissions (Figure 7). However, the 

10 fires also impacted background levels of atmospheric CO to some extent. The high VOC 

1 1 and CO emissions from the fires result in a reduction in average OH concentrations. This 

12 in turn reduces the rate of oxidation of CO and results in increased background CO 

13 levels. This increase has been estimated in the model by comparing the background CO 

14 with fire emissions (calculated by subtracting COf from total CO concentrations in BBsrf) 

15 to the CO field without fire emissions. For the Northern Hemisphere we calculate an 

16 increase in the burden of background CO of up to 1 Tg CO. 

17 5.4 Impact of Alaska/Canada Wildfires on the O 3 Budget 

18 Finally, using the model results, we examine how increases in the CO and O3 

19 concentration fields related to the wildfires affected the Northern Hemispheric and the 

20 regional trace gas budgets. In Figure 9 we show the time series for the modeled Northern 

21 Hemispheric CO and O3 burden (surface - 300 hPa) and the changes related to the 

22 emissions from the fires. The change in the CO burden reaches close to 10% around the 

23 end of July. The corresponding changes in the O3 burden are on the order of up to 4%. On 
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1 


average, the Northern Hemispheric CO and O3 burden during the summer 2004 were 


2 increased due to emissions by the fires by 4-5% and 2%, respectively. 

3 As expected, the largest changes in the atmospheric burden occurred over Alaska and 

4 Canada (50-70N, 180E-60W) with an average increase in the O3 burden of 7-9% for the 

5 altitude range surface-300 hPa, and 11-12% for the range surface-800 hPa. For 

6 comparison, over the altitude range up to 300 hPa, this is slightly smaller than the 

7 estimated contribution of stratospheric O3 (11%) in the model, but for the range up to 800 

8 hPa exceeds the contribution of stratospheric O3 (3%). 

9 Due to the transport of O3 and its precursors, effects from the fires are also expected 

10 far downwind of the source location. Over Europe (35-70N, 20W-20E) we estimate a 

1 1 contribution of O3 from the fires of up to 10% around the end of July for the surface-300 

12 hPa range and up to 8% for the range surface-800 hPa. Averaged over the summer, the 

13 contributions are on the order of 3% for both altitude ranges considered. These results 

14 show that even though the fires had a rather small contribution to the large-scale 

15 hemispheric budget of 0 3 , over certain regions and altitudes, even far downwind from the 

16 source itself, the impact is significant. 
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1 6 Conclusion 


2 We have determined the amount of O 3 produced from the wildfires in Alaska and 

3 Canada in summer of 2004 by using a combination of model simulations and 

4 observations. The modeled CO and O3 fields have been evaluated by comparison with a 

5 comprehensive set of aircraft measurements taken during the ICARTT campaign. 

6 In analyzing the O3 production from North American boreal fires we used measured 

7 and modeled CO and O3 mixing ratios at the PICO-NARE station in the Azores. The 

8 results show that the enhancement ratio AO3/ACO, defined as the increase in O3 per unit 

9 increase in CO, derived from observed and modeled concentrations at PICO-NARE is a 

10 good measure for the O 3 production from the Alaska and Canada fires. However, we also 

1 1 show that this measure can be very sensitive to the selected air masses. Our analysis 

12 yields enhancement ratios of 0.25 ppbv/ppbv for aged plumes of Alaskan and Canadian 

13 wildfires, which is in the range of values found in the literature. The enhancement ratio 

14 found for boreal biomass burning plumes is about a factor of 3-4 smaller than that of 

15 anthropogenic plumes. We have also performed a sensitivity simulation that clearly 

16 showed that the difference in the enhancement ratio for anthropogenic and boreal fire 

f7 biomass burning plumes is a result of the difference in the NO x /CO emissions ratios for 

18 these sources. 

19 Controversies exist in the understanding of the importance of O 3 production from 

20 boreal forest fires. The total net O 3 production from the boreal fires in Alaska and Canada 

21 in the summer of 2004 in our model is estimated as 6 Tg O3 which gives a contribution of 

22 about 3% to the Northern Hemispheric budget. Considering only a region spanning from 


Pfister_INTEX_revision.doc Version from 9/9/06 


Page 28/45 



1 the source to the Azores, a net chemical production of 9-1 1 Tg O3 is calculated in the 

2 model. This is in agreement with the estimate derived from enhancement ratios based on 

3 model simulations and observations at PICO-NARE (10.7-15 Tg O3). Large increases in 

4 the O3 burden are observed downwind of the fires due to transport of O3 produced near 

5 the fires as well as due to continuing O3 production in the fire plumes. Modeling studies 

6 show that the increase in the atmospheric burden of O3 is a combination of a strong O3 

7 production due to precursors emitted by the fires, and an increased destruction of O3 

8 background levels resulting from reduced peroxy radical concentrations. 

9 While the availability of satellite measurements of tropospheric CO concentrations 

10 puts constraints on the CO emissions, uncertainties remain in how to constrain the 

11 emissions of NO x and VOCs resulting in uncertainties in the estimated O3 production. 

12 Another unknown is the emissions injection height of the fires.. 

13 The results of this study indicate that fires in the boreal region can have a significant 

14 impact on the O3 production over large parts of the Northern Hemisphere. We focused 

15 our investigations on the wildfires in Canada and Alaska from summer 2004, that have 

16 been a record for this region, but comparable or even larger impacts might occur from 

17 fires in Siberia. With climate change and the possibility of increased fire activity in the 

18 Northern latitudes as a result of more frequent and/or more severe droughts and increased 

19 direct human impact [Mollicone et ah, 2006], O3 production from boreal fires might gain 

20 in importance in the future. 
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1 


Aircraft 

Species 

Instrument 

Reference 

NASA-DC8 

CO 

Tunable Diode Laser 
Absorption 

Sachse et al., 1987 


0 3 

NO Chemiluminescence 

Avery et al., submitted 


0 3 

Airborne Differential 
Absorption Lidar (DIAL) 

Browell et al., 2003 

N0AA-P3 

CO 

VUV CO Fluorescence 

Holloway et al., 2000 


0 3 

NO Chemiluminescence 

Ryerson et al. (1998) 

UK BAE146 

CO 

VUV Resonance 
Fluorescence 

Gerbig et al., 1999 


0 3 

UV Absorption 

Thermo Electron Co. Model 49 

DLR Falcon 

CO 

VUV Fluorescence 

Gerbig et al., 1996 


0 3 

UV Absorption 

Schlager et al., 1997 

MOZAIC 

CO 

Improved IR Correlation 

Nedelec et al., 2003 


0 3 

UV Absorption 

Marenco et al., 1998 


2 

3 Table 1: List of aircraft measurements included in the evaluation of the model 

4 simulations. The NOAA-P3 03 instrument is a new installation, but similar to the one 

5 described in the listed reference. 
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1 



BBsrf 

BBvert 

noBB 

BBsrf to 
noBB 

BBvert to 
noBB 

BBsrf to 
BBvert 


CO (ppbv) Mean and Std. Dev. 

CO T-test Significance 

DC -8 

1 19±37 

1 17±35 

106±31 

< 0.01 

< 0.01 

0.01 

P-3 

140±41 

1 39±4 1 

124±39 

< 0.01 

< 0.01 

0.07 

BAE146 

100±24 

99±23 

84± 1 3 

< 0.01 

< 0.01 

0.14 

Falcon 

1 10±24 

108±23 

85±1 1 

< 0.01 

< 0.01 

0.06 

MOZAIC 

1 17±39 

1 16±39 

105±36 

< 0.01 

< 0.01 

0.06 


O 3 (ppbv) Mean and Std. Dev. 

O 3 T-test Significance 

DC -8 

59± 14 

60±14 

57± 14 

< 0.01 

< 0.01 

0.02 

P-3 

59± 1 0 

59± 1 0 

57± 1 0 

< 0.01 

< 0.01 

< 0.01 

BAE146 

5 1 ± 14 

52± 1 5 

48± 12 

< 0.01 

< 0.01 

< 0.01 

Falcon 

58±1 1 

59±1 1 

53± 1 0 

< 0.01 

< 0.01 

0.01 

MOZAIC 

62±13 

63±13 

60±13 

< 0.01 

< 0.01 

< 0.01 

DIAL 

54±8 

56±9 

53±8 

< 0.01 

< 0.01 

< 0.01 


2 

3 Table 2: Statistics over the altitude range 0-8 km for the model simulations BBsrf, 

4 BBvert and noBB and the platforms listed in Table 1. Mean and standard deviation for 

5 CO and O 3 concentrations and the significance level of the Student’s T-Statistics 

6 comparing BBsrf to noBB, BBvert to 110 BB, and BBsrf to BBvert are shown. 
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n 
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0.41 
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2036. 
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— 1 3 ±3 1 

0.13 
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- 1 ±26 
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0±26 
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OX 

n = 

1802. 
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-22±26 
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-9 ±24 

0.54 
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0.56 
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1 — 14±44 
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0.61 
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0.61 
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Figure 1: Model evaluation with aircraft data for CO (blue diamonds: observations, 
black crosses: model with fire emissions injected at surface; green squares: model 
simulations with fire emissions injected over 0-9 km; red triangles: no fire emissions). 
The mean percent bias and standard deviation (model minus measurement), correlation 
coefficient r and number of data points for 2-km wide altitude bins are specified. 
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3±1 1 0.58 
n = 772. 

I.i. 1 .i.i.i.Li. 


20 40 60 80 100120 
Ozone VMR (ppb) 


DLR Falcon 


cm -11 ±16 0.58 
I~m -4*20 0.16 

rX= -5±18 0.35 
1 v 1 n=256. 

m -20±16 0.07 
cm — 14±18 -0.1 

cm — 1 5± 1 8 -0.0 

1 <> 1 n=255. 

cm — 14± 14 -0.0 
cm — 7±14 0.16 

cm -9±13 0.14 
■— V — 1 n = 192. 

-5±10 0.26 
-1 ±9 0.36 

-1 ±9 0.34 
. n=338. 


cm 

cm 

cm 

cxa 


.-L-l 


i— Lj 


20 40 60 80 100120 
Ozone VMR (ppb) 


<u 4 


M0ZAIC 


rZT~l-4±15 0.44 
con-1 ±15 0.43 

CX=~2±15 0.44 
1 V 1 n=239 1. 
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Figure 2: As Figure 1, but for O 3 . 
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Figure 3: Measured and modeled time series of CO and O 3 mixing ratios at PICO- 


NARE (daily averages are shown). 
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Figure 4: Observed (a) and modeled (b,c) CO-O 3 relationships at PICO-NARE. Air 


masses are separated into three groups: mostly non-fire related origin (blue), some 


biomass burning impact (green), pronounced biomass burning impact (red). In (c) results 


are also shown for a simulation where the NO x fire emissions were increased by a factor 


of 10. 
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Figure 5: Observed (a) and modeled (b,c) CO-O 3 excess mixing ratios. In (a) and (b) a 
constant background (bg) is calculated from the subset of non-fire impacted airmasses, in 
(c) the background is derived from the noBB simulation. 
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2 Figure 6: Enhancement Ratio (mean slope and standard deviation) determined by the 

3 scatter and the enhancement technique as a function of intensity of biomass burning 

4 influence of considered air masses. The number of selected data points is represented by 

5 the shaded area. Results for the BBsrf simulation. 
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3 Figure 7: Correlation between the model fire tracers COf and 0 3 NOx and the difference 

4 in CO and O 3 mixing ratios between the simulations BBsrf and noBB (defined as dCO 

5 and d0 3 ) at PICO-NARE. 
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1 Figure 8: Statistics for ozone concentrations and ozone production and loss terms for 

2 three different subsets of fire plume intensity as characterized by the ratio of COf/CO. (a) 

3 Spatial distribution for selected data points, (b) Frequency distribution for the height of 

4 the selected data points (c) Frequency distribution of the volume mixing ratios O 3110 BB 

5 and O3. (d) Net O3 production rate for the fire tracer 0? U,A (e) Frequency distribution of 

6 the 03 production rate for O^noBB and OS (background O 3 estimated by subtracting O3 

7 and 0/ VOa ). (f) as (e) but for the chemical lifetime (g). Frequency distribution of the 

8 volume mixing ratios 03 NOx , (IO3 and OS. 
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3 Figure 9: Northern Hemispheric burden of CO and O 3 for the altitude range surface-300 


4 mbar. Model results for simulations BBsrf (solid line) and noBB (dotted line). The dash- 


5 dotted line denotes the percentage difference between simulations BB and noBB (shown 


6 on the secondary ordinate). 
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